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Abstract 

We examine in detail the predictions of the charged majoron model, introduced re- 
cently by Burgess and Cline, for + — > + j3j3 transitions. The relevant nuclear ma- 
trix elements are evaluated, within the quasiparticle random phase approximation, for 
76 Ge, 82 Se, 100 Mo, 128 Te and 150 Nd nuclei. The calculated transition rates turn out to be 
much smaller than the experimental upper limits on possible majoron emission, except in 
a small region of the model's parameter space. 
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In 1987 Elliott, Hahn and Moe JI| observed, for the first time, the two-neutrino double 
beta decay (Pfou) of 82 Se into 82 Kr by a direct counting method. Almost simultaneously, 
Vogel and Zirnbauer showed that, within the quasiparticle random phase approxima- 
tion (QRPA), it was possible to explain the smallness of the measured Pfov decay rates. 
Since then, impressive progress has been achieved in the experimental investigation, and 
the 2i/-decay mode has been unambiguously observed in several nuclear isotopes ||. This 
process occurs at second-order in the charged-current weak interactions in the Standard 
Model, and is the slowest process measured so far in nature. As such, it offers a unique 
opportunity to test the nuclear structure models for half-lives ^ 10 20 years. 

However, the renaissance of interest in /?/?-decay over the last decade is mostly stim- 
ulated by the possibility of observing other decays to which these experiments are also 
sensitive. The hope is to find a smoking gun for "new physics" from beyond the stan- 
dard model. The most promising processes of this type are the lepton-number violating 
neutrinoless decay ((3{3q v ), and the decay (3/3m, in which the two outgoing electrons are 
accompanied by a Nambu-Golstone boson, called the majoron. Both processes were pre- 
dicted ||, |5|] by the model introduced by Gelmini and Roncadelli ||. While this simple 
and elegant model stimulated many experimental searches, it was subsequently found to 
be incompatible with the LEP measurement of the invisible width of the Z boson |7], H|. 

Thus, to account for an anomalous excess of /3/3 events for which some experimental 
groups had preliminary evidence ||, Burgess and Cline recently advocated a new class 
of "charged majoron" models [ID, ll| so-named because their majoron carries the U(l) 
charge of lepton number, presumed to be unbroken. Refs. [10, |TTJ estimated the nuclear 
matrix elements needed for charged majoron emission to account for the anomalous events 
observed in 76 Ge, 82 Se, 100 Mo and 150 Nd nuclei. In the present work, we compute the 



transition rate in detail for the particular model given in [pJ]] , simplifying their analytic 
expressions, and evaluating the corresponding nuclear matrix elements for the above- 
named elements as well as for 128 Te. We use the quasiparticle random phase approximation 
(QRPA) [|l|, [131 ) which has been shown to give good estimates for the flfyv transition 
probabilities. Finally, the resulting transition probabilities are compared with the present 
experimental data. 



In the charged majoron model (CMM) of ref. [jO], the standard model gauge symmetry 
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group is augmented by a global nonabelian flavour symmetry group, SUf(2) x Ul>(1), 
which breaks down to the ordinary lepton number Ul(1) subgroup (see also refs. [|14], |T5||). 
To implement this symmetry-breaking pattern, an electroweak-singlet, 5'L r F(2)-doublet 
scalar field $ is introduced, which gets a vacuum expectation value. The model also 
includes nonstandard electroweak-singlet Majorana neutrinos (N + ,N_), s + and s_, and 
the resulting lagrangian density which respects all the symmetries is 

C = -\LHP r s_ - Ms + P R s_ - g + (NP L s + ) $ - g. (iVP L s_) <£ + h.c. (1) 

Here Pl and Pr denote the projections onto left- and right-handed spinors, L is the usual 
lepton doublet, and H is the charge conjugate of the Higgs doublet. After symmetry 
breaking by the vacuum expectation values (H) = v = 174 GeV, ($) = u ~ 100 MeV, 
the resulting neutrino mass matrix yields a massless neutrino v' e and two heavy Dirac 
neutrinos ip± with masses 

1 r. 



.U- M 2 : \JM^ - g 2 + u 2 (\ 2 v 2 + g 2 u 2 ); M = - \M + X v + (g 2 + + gi)u A \ . (2) 
In terms of the mass eigenstates, the electroweak eigenstate is 

We) = C d \v' e ) + S d C a \lp_) + S 9 S a \i) + ) (3) 

where sg = sin 8, cq = cos 8, etc., denote the mixing angles with 



\ v 2MJ\ 2 v 2 + g 2 u 2 

tan# = ; tan2a = — v r^. (4) 

g_u M 2 — \ z v z + (gj_ — gi)u z 

The transition probability for the + — > + /3/3-decay has the form [JTT ] 

T(m = (G 256 S J C)4 lAm |2 S {Q ~ 61 " 62)3 A fo^fe)** (5) 



in the notation of ref. [O], with the transition amplitude given by 

A{(3(3) = 2V2(s 2a S28S8)ig + g 2 A (0 + \\^2Y Rnm ■ r nm dh a /dr nm \\0 + ). (6) 



Here 



(27r) 4 (fcg - /r + ze) 
P t = (e-M' + te)- 1 ; M = 0, (7) 
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with r mn and /i being, respectively, the separation in the position between two decaying 
nucleons and the average excitation energy of the intermediate nuclear states, and 

1 2 



Y 



Rnm 



:(T n C m — C n er m ) + i 



#1 



9a 



(<x n x D m + D n x a m ) + 



9y_ 

9a 



;D n - D m ) 



happens to be same nucleon recoil operator as appears in Doi et al, |J, except for the 
sign of the second term. (We believe the sign difference is due to an error in ref. 

It is possible to simplify the above expressions for the matrix element. Previously 
it has been argued that the middle term of eq. (§) is the most important. Within the 
approximation where the momenta of electrons and majoron are neglected in comparison 
with those of the nucleon, and because the main contribution involves the spin singlet 
state of two nucleons, this term gives || [l 



* Rnm ' f"r 



2M, 



TV 



(<T n " Cm) (^ntn ' f° nm ), 



(9) 



where M N is the nucleon mass, fa ~ 5.6 (for g A = 1 [|T^, and V nm = d/dr nm . 
From eq. ([]) it is clear that the amplitude is largest if S2 a = 1, so we shall make that 
assumption. Furthermore the mixing angle 9 is typically constrained to be small, so we 
take sg = 9. Then, after performing the po integral in eq. (|7]) and using the identity 
V ■ vdh a {r) / dr = V 2 /i Q (r), the amplitude can be written as a difference of two pieces, 



V2n 



M + - Mr 

J l CM J l CMl 



where 



[r 

2M, 



N 



(0 + \\J2h(r mn ; M±) <r n ■ cr m ||0 + ), 



(10) 



is the nuclear matrix element for charged majoron emission corresponding to the exchange 



of the neutrino with mass M±, and 



h(r mn ; M ± ) 



Ml 



[h(r mn ; M ± ) - h(r mn ; 0)] + 



1 d 



2d Ml 



h(r. 



is the corresponding neutrino potential, with 



h(r mn ; M ± ) 



dk 

2^ 



k 2 



UJ 4 



k 2 + M_ 



1/2 



(12) 



(13) 
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The inverse half-life can be now cast in the form 



[ T (o+ -> 0+)]- 1 = g 2 CM g CM \M CM \\ (14) 

where 

9cm = g + e 2 /2, (15) 
is the effective majoron coupling, and 

(G F g A cos Q c ) 



1287r 7 hi2 



f(Q-e 1 -e 2 ) 3 f[k i e i F(e i )de i , (16) 

i=i 

is the kinematical factor as defined in ref. ||. 

To perform the nuclear structure calculation we use the Fourier-Bessel expansion of 
the charged majoron matrix element. Thus 

Mt M = E m(M ± ;L,S,J*), 

LSJ-* 

where L, S, J and tc are, respectively, the orbital angular momentum, the spin, the total 
angular momentum and the parity of the intermediate nuclear states. Within the QRPA 
formulation presented in ref. [T2[ the individual nuclear matrix element is 



m(M ± ; L, S, J n ) = (-) s £ W™ J W$ r Tl L (pnj> , n'; M ± )A + (pn; aJ^A* (p'n'; a J*). 

apnp'n' 

(17) 

The amplitudes A ± (pn;aJ 7r ) are 

A+(pn; aJ) = ^p p p Q [u p v n X pn . aJ + v p u n Y pn . aJ ] ; 
A_(pn; at J) = ^p p p n [v p u n X pn . aJ + u p v n Y pn . aJ } , 

where the unbarred (barred) quantities indicate that the quasiparticles are defined with 
respect to the initial (final) nucleus; p p x = u 2 + v 2 p^ 1 = u 2 l + v 2 , and all the remaining 
notation has the standard meaning |T2[] . The angular momentum and radial pieces in flTT| ) 
are, respectively, 

1 q 

-n 9 jn 



W^ n SJ = ^ +L V2JSLj p ]J n (£ n 0L0\£ p 0){ L S J 

'P 2 ^P 
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and 

roo 

ft L (pnp'n';M ± ) = J dkk 2 v(k;M ± )R^ n (k)R^ Ql (k), 



with 



9 k 2 

v(k;M ± ) = — 



IT [Ml 



uj±{lo± + fi) k(k + fi) 



1 d 
+ -x 



2 3M|w ± (a; ± + /i) 



R L( k ) = / u n (r)u p {r)j L (kr)r 2 dr, 



and 

it(r) being the single-particle radial wave functions. 

The numerical calculations were performed with the 5-force (in units of MeV fm 3 ) 
V = —4tt(v s P s + v t P t )5(r), with different strength constants v s and v t for the particle- 
hole, particle-particle and pairing channels. An eleven-dimensional model space was used, 
including all the single particle orbitals of oscillator shells 3hu and AHlj, plus the O/2.9/2 and 
O/tn/2 orbitals from the 5hu> oscillator shell. The single particle energies, as well as the 



parameters vf atr (p) and vf a "*(n), have been fixed by the procedure employed in ref. [12 
(i.e., by fitting the experimental pairing gaps to a Wood-Saxon potential well). 

The dependence of M.% M on M ± is illustrated in Fig. 1 for the 76 Ge— s> 76 Se decay. 
As one might expect, the matrix element is insensitive to the mass until it starts to 
exceed the Fermi momentum of the nucleons, around 100 MeV, thereafter giving a 1/M 2 
suppression. Since the matrix element M. C m is the difference between the M + and M_ 
contributions, one can easily recover the result for an arbitrary pair of masses by taking the 
difference between the two corresponding matrix elements. Obviously for M_ = M + there 
is destructive interference between the matrix elements M-cm arLC ^ M.% M . This happens, 
for example, with the choice of parameters M = Xv and g± ~ 1, which from eqs. (0) and 
(||) yields M ± = g + u\/\ ± 9. Since the mixing angle 9 is experimentally constrained to be 
of the order of 0.1, [] the calculated values of [T^^] -1 turn out to be four to five orders 
of magnitude smaller than the corresponding experimental upper limits for the majoron 
emission, due to the additional 9 2 suppression. These limits are displayed in Table 1 for 
the 76 Ge, 82 Se, 100 Mo, 128 Te and 150 Nd nuclei. For the sake of completeness, the table 
also shows the measurements for the ftfov decays, and Q CM for the effective axial-vector 



coupling g A — 1 |T2|, |17 |. 

1 For a discussion of these experimental constraints, see ref. p0| . 
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Figure 1: Charged majoron matrix element M.% M (in natural units) for the 76 Ge— > 76 Se 
decay, as a function of the heavy neutrino mass M ± (in units of MeV). 

Figure 1 also shows that the most favorable situation for majoron emission occurs for 
in limit that the heavier neutrino becomes infinitely massive, thus making no contribution 
to the total rate. For M_ we will use 100 MeV, since this is the largest value which still 
gives an unsuppressed amplitude. In Table 1 we compare the theory to the data by 
showing how large the effective coupling g CM would have to be in this case in order for 
the CMM rate to be equal to the experimental limit on the majoron-emitting mode of 
(3(3 decay. For all the elements, g CM must be 0.1 — 0.2 to be presently observable. Given 
the above mentioned experimental constraints 9, from eq. fllED it is clear that one would 
need a strong coupling in the sterile neutrino sector in order to achieve such a large value 
of the majoron-emitting rate. 
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Table 1: The size of the effective coupling g CM which would be needed for the rate of 
emission of charged majorons to be equal to the experimental limit on this process; Q CM 
and |.M| CM are the corresponding kinematical factors and the nuclear matrix elements, 
respectively. For the sake of comparison the pertinent experimental data for the two- 
neutrino processes are also shown. 



Nucleus 






GcAi(yr~ 


- 1 ) 


\M\ CM 


g C M needed 


7e Ge 


(6.99 ±0.08) 10~ 22a ) 


< 6.0 10" 23 a ) 


2.10 10" 


-20 


0.28 


< 0.19 


82 Se 


(9.26+°-|) 10- 21fe ) 


< 6.2 10~ 22 c ) 


3.55 10" 


-19 


0.28 


< 0.15 


W0 Mo 


(8.69^;p) 10- 20rf ) 


< 3.0 10~ 21 e ) 


7.33 10- 


-19 


0.31 


< 0.21 


128 Te 


(1.30 ± 0.07) 10 - 25/ ) 


< 1.3 10- 25 ^) 


5.24 10- 


-22 


0.27 


< 0.059 


150 Nd 


(5.88±1;1^ ± 1.21) KT 20 ^ 


< 1.9 10" 21 ^ 


6.06 10" 


-18 


0.13 


< 0.077 



a ) (laboratory data) ref. [18 1 

b ) (laboratory data) ref. |19|] 

c ) (laboratory data) ref. J]J 

d ) (laboratory data) ref. p0| 

e ) (laboratory data) ref. [^] 
*) (geochemical data) ref. [22 
9 ) (laboratory data) ref. P3 
h ) (laboratory data) ref. [|4]] 



In conclusion, we have found that the rate of majoron-emitting, neutrinoless (5(5 decay 
in the charged majoron model is unobservably small unless there is large mixing between 
exotic sterile neutrinos of mass ^ 100 MeV, and strong couplings among the sterile 
neutrinos. In computing the relevant nuclear matrix elements, we have not considered 
the effects of finite nucleon size or short-range two-nucleon correlations, which would tend 
to reduce the calculated matrix elements |12|]. On the other hand, the arguments used 
to simplify the nuclear transition operator from (Q) to @ are not rigorous, and it is also 
possible that future variants of the CMM considered here might evade the suppression of 
the rate by the mixing angle 9. Thus, while we believe that ours is the most quantitative 
analysis of the CMM to date, if the experimental situation should give serious indications 
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of anomalous (3(3 decay events in the future, it would become appropriate to undertake a 
yet more careful evaluation of the model's predictions. 



Note added: As we were finishing this work, Hirsch et al. presented results in- 
cluding a somewhat less detailed analysis of the CMM, in which they reached conclusions 
similar to ours. 
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